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Background It has been recently shown, that some Skyrme functionals can lead to non-converging results in the calculation 
of some properties of atomic nuclei. A previous study has pointed out a possible link between these convergence problems 
and the appearance of finite-size instabilities in symmetric nuclear matter (SNM) around saturation density. 

Purpose We show that the finite-size instabilities not only affect the ground state properties of atomic nuclei, but they can 
also influence the calculations of vibrational excited states in finite nuclei. 

Method We perform systematic fully-self consistent Random Phase Approximation (RPA) calculations in spherical doubly- 
magic nuclei. We employ several Skyrme functionals and vary the isoscalar and isovector coupling constants of the 
time-odd term s • As . We determine critical values of these coupling constants beyond which the RPA calculations do 
not converge because RPA the stability matrix becomes non-positive. 

Results By comparing the RPA calculations of atomic nuclei with those performed for SNM we establish a correspondence 
between the critical densities in the inhnite system and the critical coupling constants for which the RPA calculations 
do not converge. 

Conclusions We hnd a quantitative stability criterion to detect hnite-size instabilities related to the spin s • As term of a 
functional. This criterion could be easily implemented into the standard htting protocols to hx the coupling constants 
of the Skyrme functional. 

PACS numbers: 26.60.Gj, 21.60.Jz, 21.30.-x 


I. INTRODUCTION 

The Nuclear Energy Density Functional (NEDF) is a 
standard tool used to describe properties of atomic nu¬ 
clei from drip-line to drip-line [l|, . Among the differ¬ 

ent functionals available nowadays, the one derived from 
the non-relativistic Skyrme interaction is one of the most 
popular. A very important aspect for the determination 
of a functional is the optimization procedure used to ob¬ 
tain its effective coupling constants il. To this pur¬ 
pose one has to minimize in a multi-dimensional space a 
penalty function. The latter is obtained by considering 
some observables of finite nuclei as for example masses 
and radii and some pseudo-observables of infinite nuclear 
matter [^-Q. The minimization procedure can be quite 
evolved since the penalty function presents discontinu¬ 
ities in the parameter space and therefore may induce 
some difficulties in finding a stable minimum. For a more 
detailed discussion, we refer the reader to Ref. @ and ref¬ 
erences therein. 

Since the input parameters may present some uncer¬ 
tainties or their number is limited, a coupling constant 
can be poorly constrained making the extrapolation pro¬ 
cedure very problematic M- Such a situation can also 
eventually lead to unphysical instabilities. By instance, 
in Ref. [^, Margueron et al. have shown that most 


of the Skyrme functionals present spin or ferromagnetic 
instabilities, that is a spontaneous and complete polar¬ 
ization of the infinite medium. The latter is not pre¬ 
dicted by any other model either phenomenological or 
ab-initio thus leading to the conclusion that it is a 
pathology of the Skyrme functional itself. These type of 
instabilities can be simply controlled by inspecting the 
properties of the corresponding Landau parameters de¬ 
rived from the Skyrme functional [l3l| during the opti¬ 
mization procedure [l^ . However, another type of in¬ 
stability, not related to the entire system, but specific to 
finite-size domains with typical size A ss 27 r/( 7 , where q is 
the transferred momentum has been recentl y de tected by 
Lesinski et al. [l^ for the SKP [T^ and LNSfr^ function¬ 
als. Contrary to the previous case, since the instability 
occurs at non-zero values of q, it can not be detected by 
exploring the Landau limit. 

In Ref. [1^, Hellemans et al. investigated the nature 
of finite-size instabilities for the scalar-isovector channel. 
The authors have established a simple relation between 
the presence of these instabilities in finite nuclei calcu¬ 
lations and the position of zero-energy modes calculated 
using Random-Phase-Approximation in infinite nuclear 
matter M- In particular, it has been shown that when¬ 
ever the critical density pc at which a pole {i.e zero en¬ 
ergy mode with infinite strength) appears in SNM is close 
to the saturation density of the system, that instability 
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manifests itself in finite nuclei leading to non-physical so¬ 
lutions. For example, in the scalar-isoscalar case it leads 
to a complete separation between protons and neutrons. 
Such a relation is particularly important since RPA cal¬ 
culations in SNM [I^ are analytical and not time con¬ 
suming. Recently, they implemented a criterion based 
on RPA in SNM into a fitting procedure to prevent, by 
construction, this kind of problems [^ . 

In the present article, we explore deeper the question of 
the correspondence between poles in the infinite medium 
and instabilities in finite nuclei. In particular, we want 
to know whether this remains valid in the vector channel 
of the interaction. Our study is motivated by the prob¬ 
lems detected recently using Time-Dependent-Hartree- 
Fock methods (TDHF) [2l| and Cranked-Hartree-Fock- 
Bogoliubov (CHFB) calculations [1^, which clearly ex¬ 
hibit problems of convergence when using some Skyrme 
functionals. In particular, in Ref. [l^, the instability 
problem has been clearly associated with the time odd 
terms which depend on the laplacian and gradient of 
the spin density [2^. In the present article, we follow 
the protocol outlined in Ref. [l^, but using the new 
spherical RPA code recently developed by the FIDIPRO 
group [i^, [13 • The main advantages of using RPA are 
that the time-odd terms [13 are active and that it is 
numerically less time-consuming and thus more adapted 
for systematic calculations as compared to other methods 
such as CHFB or TDHF. 

The article is organized as follows: in SeclTT]we briefly 
sketch the formalism of the Skyrme functionals, in Sec IHIl 
we present the method used to detect and identify the 
appearance of finite-size instabilities in finite nuclei, while 
in Sec lIVI we present the formalism of RPA in infinite 
matter. The link between the two methods is presented 
in Secl3 where a simple criterion to detect instabilities 
is also proposed. Finally we present our conclusions in 

SecED 


II. SKYRME FUNCTIONALS 

The Skyrme functional, £sk, is a linear combination of 
coupling constants and local densities of the form 

£sk = + Ct*[po]st+Cf^PtApt 

+ ■ Sf)^ -|- C(^®stAst -|- Cl [ptTt — jt ) 


and can be found in Refs.[^.[^. The density dependent 
coupling constants are written as 


CM = Ct^ + p<^cr, (2) 

CM = Cf + pZCr. (3) 


The finite-size instabilities are mainly related to the 
presence of derivative terms: pt^pt for the scalar chan¬ 
nel and StAst for the vector one. The term C^^C^ ■ St)^ 
represents also a possible source of instability as already 
mentioned in Ref. [22|, but it is active only if the inter¬ 
action contains an explicit tensor term. In the present 
analysis, we do not consider explicitly the instabilities 
related to such a term and we focus more on the ones 
related to the time-odd Laplacian term. 

To perform our study, we chose six representative 
Skyrme functionals: T44 [23L SAMi [il], BSk27 [^ . 
SHI [ 33 , SKO' [m and SLybQ . These functionals have 
been obtained by different groups using different penalty 
functions and they are currently used within the nuclear 
structure community to perform calculations of nuclear 
observables for both ground and excited states of nuclei 
along the nuclear chart. 

The presence of long-wavelength ferromagnetic insta¬ 
bilities in SNM can be easily determined by inspect¬ 
ing the behavior of the Landau parameters C?^=o,i and 
Cf^Q 1 , which can be expressed in terms of Skyrme cou¬ 
pling constants as [13 


No^Co = 2C*’° + (2 + a)(l + a)Co^’“p^ 

+ 2kM + > (4) 

Aq-IGi = -2klCl - , (5) 

A-^G; = 2Cf + {2 + a){l + a)CrPo 

+ 2kM + \klCf , (6) 

Aq-'g; = -2klcl - “^klCf , (7) 

where kp the Fermi momentum and Nq = is 

the density of states at the Fermi surface. In Figlll we 
show the evolution of the Landau parameters G; and GJ 
in SNM for the six selected functionals. 

According to the results of Ref. m , whenever the in¬ 
equality 


Cj s, . T* - 


JtpuJtpu 


Gi>-{21^1), ( 8 ) 





+ G7''(ptVJi+stVxjt)} . 


The definition of the local densities pt,St,... is standard 


is not satished, a static deformation of the Fermi surface 
against spontaneous polarization is observed. In pres¬ 
ence of an explicit tensor term new stability condition are 
necessary (see Refs. [H, [13 for a more detailed discus- 
(I(gion). From Fig.[l] we observe that the SKO' functional 
is spin-unstable in the isoscalar channel at very low den¬ 
sities p Ri 0.14 fm“^. Such a problem manifests in RPA 
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TABLE I: Values of the coupling constants C^' 0,1 and SNM 
saturation density p^at of the functionals considered here. 



FIG. 1: (Color online) Evolution of Landau parameters Gi 
and G'l {I — 0,1) as a function of the density in SNM for the 
six Skyrme functionals considered in the article. 


calculations of finite nuclei. We have selected the spher¬ 
ical doubly-magic nuclei ^°Ca, ®®Ni, ^^^Sn, and as 

specific examples. In particular we have not been able to 
perform any complete calculation due to the appearance 
of several imaginary phonons which spoil the results. It is 
actually possible to obtain a result only by reducing the 
available ph-space, but in that case the results strongly 
depend on the parameter space and on the related cut¬ 
off. The use of small parameter space is the reason why 
other groups have been able to present results using this 
functional, although the physical value of such results is 
very doubtful. For such a reason the SKO' should not be 
used to perform RPA calculations. 

The functionals SLy5 and T44 are also unstable, in the 
spin-isovector channel, but at much larger values of the 
density which are not relevant for finite nuclei, although 
they can be important for astrophysical applications [3^ . 

It is worth noticing that the functionals SAMi and 
BSK27 have been built to remove these instabilities di¬ 
rectly during their fitting procedure. Therefore, they re¬ 
spect Eq. ([5]) up to several times saturation density. 

The functionals Sill and BSk27 give Landau parame¬ 
ters with I = 1 identically equal to zero. This is due to 
the dropping of the so-called term in these function¬ 
als, i.e. we have = 0. By inspecting Eqs. (Himi . 

we observe that Landau parameters do not depend on 
the coupling constants and As discussed in 

Ref. , these terms are proportional to in the resid¬ 
ual interaction, and thus they are zero in the limit g —> 0 . 


SIIII 

SLy5 

T44 

BSk27 

SAMi 

SKO' 


Psat(fm-^) C^®(MeVfm'=) Cf®(MeVfm^) 


0.145 

0.160 

0.160 

0.158 

0.159 

0.160 


17.031 

46.087 

72.670 

27.925 

44.720 

42.790 


17.031 

14.062 

1.393 

16.691 

20.962 

16.553 


III. RPA INSTABILITIES IN FINITE NUCLEI 

We now focus our attention on the S(As( terms, to 
analyze its effects on possible instabilities. The value of 
the coupling constants is given in Tab. U for the 
functionals considered in this article. 

To detect the presence of finite-size instabilities, we 
perform RPA calculations in the selected set of nuclei 
using the numerical code described in Ref. [ 2 ^ . The 
calculations are fully self-consistent, i.e. include both all 
the terms of the functional and all the states within the 
basis. In the present article we consider a basis of 16 
major oscillator shells. No further cut-offs are included 
in the calculations as commonly done, see for example 
Refs. (35l - [^ for discussion. It is important to stress 
that although the use of a cut-off in ph space can be 
safely adopted for stable calculations, this could hide the 
problems discussed here. To observe the presence of in¬ 
stabilities, we calculate the stability matrix S [s^ 


5 = 


A B* 
B* A* 


(9) 


which should be positive-defined in the case of RPA. Neg¬ 
ative or imaginary eigenvalues of the S matrix are thus 
the signal of finite-size instabilities, which are related to 
the fact that the ground state configuration we have used 
to calculate the vibrational band is not a minimum of the 
energy surface. As a consequence, the system is no more 
stable against small fluctuations around this point. To 
study the impact of the term SqAso, we put to zero the 
coupling constant and we vary C^® and vice versa 
through a multiplicative parameter 7 . Similarly to what 
has been done in Ref. [l^ for the term piApi, we want 
to determine the critical value of the coupling constant 
beyond which imaginary eigenvalues appear. 

In Fig. [2] (a), we show the evolution of the lowest 
phonon energies in ®®Ni as a function of the parameter 7 
which multiplies the isoscalar coupling constant for the 
SAMi functional. For 7 = 1 we obtain the nominal value 
of the coupling constant as shown in Tab. [H When we 
modify 7 , we are effectively building a series of different 
functionals which differ only by the value of this time- 
odd coupling constant, and for each functional we solve 
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FIG. 2: (Color online) Evolution of the lowest energy vibra¬ 
tional states for different multipolarities in ®®Ni as a func¬ 
tion of the isoscalar coupling constant 0^“ (panel a) and the 
isovector one (panel b), for for the SAMi functional. 


the corresponding RPA equations, obtaining the energies 
of excitations for multipolarities up to J = 6 and both 
parities tt = ±1. All these newly constructed functionals 
have exactly the same bulk properties of infinite matter, 
ground state properties of even-even nuclei and Landau 
parameters. 

When the multiplicator is in the range 7 £ [0 : 0.8], the 
phonon spectrum is essentially flat, telling us that the 
phonons are not so sensitive to the exact value of this 
coupling constant, but when we approach 7 « 1 (be. 
the nominal value), we notice a sudden bending in the 
J = S’*" phonon. Going beyond the nominal value, we 
observe that there is a change in the ordering of the low¬ 
est excited states until we reach a critical value at which 
the RPA calculations are no more trustful since some of 
the eigenvalues become imaginary. The critical value is 
not the same for all different multipolarities, still such a 
value exists for all multipolarities and we observe a con¬ 
centration of the strength of the low-lying RPA phonons 
when the respective critical value is approached. We re¬ 
mind that the breathing mode 0 + does not break the 
time-reversal symm etry thus the time-odd coupling con¬ 
stants play no-role (^ . 

We have also analyzed the behavior of the vibrational 
states up to J = 12 and observed the same trend. It also 
seems that they are not decisive for setting the critical 
value of the coupling constant. For this particu¬ 
lar parameterization, the isoscalar coupling constant is 
closer to its critical value for most of the multipolarities 
considered here. 

In Fig. [5] (b), we repeat the same calculations, but for 
the isovector coupling constant The multiplier pa¬ 

rameter 7 is set in exactly the same way, but for the 
isovector constant. In this case the critical coupling con¬ 
stant is located at Ri 2.5 its nominal value. 


In Fig.[3l we investigate further the instabilities related 
to the isoscalar coupling constant in ®®Ni, by perform¬ 
ing systematic RPA calculations with the other selected 
functionals. We observe that T44 is not stable when we 
use the nominal value of the coupling constant, making 
it improper for RPA calculations. To converge the RPA 
calculations, we have to use a smaller multiplier 7 < 0.5. 
Since in the present study, we do not modify the cou¬ 
pling constants related to the term (Vst=o,i)^j our result 
can not be directly compared with the one obtained in 
Ref. [Il|. 

Similarly the isoscalar coupling constants of SLy5 and 
BSk27 parameterizations of the Skyrme force are close 
to the limit of stability and we clearly observe that the 
phonons in the area where 7 r; 1 strongly depend on the 
exact value of the time-odd coupling constant We 

can thus conclude that these functionals are not adapted 
to describe vibrational states in finite nuclei. 

We have also tested the dependence of our results on 
the size of the basis. In Fig. [4l we show the evolution 
of the lowest critical coupling constant (panel a) 

and G^® (panel b), i.e. the value beyond which at least 
one of the multipolarities studied here gives an imaginary 
phonon, as a function of the number of major shells n 
included in our calculation for the Sill functional in ®®Ni. 
Since we perform the calculations on a discrete mesh of 
the multiplicative factor 7 , we have defined the point of 
collapse as the value in between the last point ^n-i for 
which all the phonon energies are real and the first point 
7 Ar at which we have a collapsed solution. We thus define 
7c = [iN + 7 Af-i )/2 with the error bars set by the values 
of 7 Ar and 7 Ar_i. 

We observe a strong dependence of the coupling con¬ 
stant G^® on the number of shells used in the calcu¬ 
lations: the larger n, the smaller G^®. Such a feature 
has been also observed for the coupling constant in 
Ref. [ 1 ^. From this figure we can argue that the use of 
very small basis or extra cut-offs on the ph-space in RPA 
calculations could hide this problem. 

It is hard to judge what extrapolation one could make 
to properly obtain the critical value of the coupling con¬ 
stant at n ^ 00 . For such a reason, we performed an 
extrapolation by using three functions of the form 


/.(n) = ^+5,, ( 10 ) 

with i = 1,2, 3. Using the proposed functions, we have 
set the lower limit of the coupling constant by fitting the 
function fi and the upper limit fitting the function /a on 
the calculated data. The fit done with the function /2 
obtains the most reliable critical value, as this function 
gives the best value when fitted to the calculated data 
points for all of the parameterizations of the Skyrme force 
used in this article. 

Another important aspect is the dependence of the po¬ 
sition of the critical coupling constant for the considered 
nuclei. It is important to notice that as in Ref. [l^ , shell 
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FIG. 3: (Color online). Evolution phonons in ®®Ni as a function of the multiplicative factor 7 for T44 (panel a), SLy5 (panel 
b), BSk27 (panel c) and Sill (panel d). The legend is the the same as for FigO 
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TABLE II: Critical parameters, obtained for some Skyrme 
parameterizations in the limit of infinite number of shells. 



^Oc 

(MeV fm)® 

r>^s 

^Ic 

(MeV fm)® 


lower 

best 

higher 

lower 

best 

higher 

Sill 

27.28 

32.37 

34.05 

36.55 

43.10 

45.25 

SLy5 

41.26 

47.69 

49.80 

36.90 

43.47 

45.63 

T44 

39.06 

47.28 

49.96 

-13.02 

-2.01 

1.60 

BSk27 

22.11 

28.22 

30.05 

33.16 

38.81 

40.67 

SAMi 

41.30 

47.54 

49.58 

38.81 

45.15 

47.23 


FIG. 4: (Color online) Evolution of the critical coupling con¬ 
stants for the Sill functional for different number of shells 
included in the calculation. With lines we have shown some 
proposed functions to describe the dependence of the critical 
coupling constants with the shell numbers. 
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FIG. 6: (Color online) Evolution of the low-lying spectrum in 
^®^Sn (calculated with 16 oscillator shells) with respect to C'o'® 
for the SLy5 functional. In panel (a) we present the spectrum 
calculated using C^"=9.22 MeV fm®, panel (b) C^'’=46.10 
MeV fm®, panel (c) C'^“=55.3 MeV fm® and panel (d) (7^“= 
59.91 MeV fm®. 


FIG. 5: (Color online) Ratio among the critical coupling con¬ 
stant and its nominal value using the best extrapolation as 
defined in the text and for the different nuclei considered in 
our analysis. 


effects could hide or enhance the appearance of such a 
kind of instabilities. In Fig. [5] we represent the ratio 
/Cq‘^ for RPA calculations using the best extrapola¬ 
tion and for the four nuclei considered here. We notice 
that in all cases the N = Z nuclei are the most sensitive 
to the appearance of these instabilities, as discussed in 
Ref. [l^, the particular shell structure could be at the 
origin of this phenomenon. 

In Table HI we summarize the results of our analysis, 
extrapolation for infinite basis size. 

In our analysis, we pay particular attention to the be¬ 
havior of the low lying dipole excitations, since recent 
results show the appearance of pigmy dipole reso¬ 
nances in atomic nuclei. Moreover low-lying dipole states 
play an important role in describing some relevant astro- 
physical process. We recall, that in our calculations, all 


the RPA phonons in the J = 1“ channel have been or- 
thonormalized against the spurious state by construction 
as explained in Refs, [l^, |4l|, thus our results are not 
polluted by spuriosity. Varying the isoscalar coupling 
constant, we have calculated the strength concentrated 
in the low-lying states in ^^^Sn using the SLy5 parame¬ 
terization and different nominal values for the isoscalar 
coupling constant In Figl5]we show the results. It is 
clearly seen that when the value of the coupling constant 
gets close to its critical value, we observe a noticeable 
shift of the strength to lower energies, and there is some 
fragmentation of this strength over several states. This 
result seems to be artificial and marks a possible weak¬ 
ness for some calculations performed for the distribution 
of the low-lying strength using functionals with coupling 
constants that are close to their critical value. It is worth 
emphasizing that once we are above the critical value, 
the energy weighted sum rule can not be fulfilled any¬ 
more since some of the strength is taken away by the 
imaginary phonons. 

Finally we stress that even though some numerical ef¬ 
fects could influence the precision of our calculations, the 
effect of collapsing solutions in a very tight region of 
values, could be reproduced using different approaches 
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realized in different numerical codes. In particular the 
results obtained with the present RPA code has been 
checked [i^ using another RPA code [i^ based on Fi¬ 
nite Amplitude Method (4ll| . 


IV. RPA INSTABILITIES IN INFINITE 
MATTER 

In this section, we present the formalism of the linear 
response theory in inhnite nuclear matter. This method 
has been the subject of a recent review article , where 
all the details of the formalism are presented and dis¬ 
cussed. We limit ourselves to sketch the basic ingredients 
of the formalism. 

The first ingredient is the Hartree-Fock retarded prop¬ 
agator of a non-interacting ph-pair Since for the present 
study we ignore charge-exchange process, the particle 
and the hole in the same pair share the same isospin 

I 


number r = p, n 


9 {k-p-k)-e{k^p-\k + q\) 




LO + eT-(k) - er(k 4- q) -I- f ?7 


( 11 ) 


where 0 {kp—k) is the standard step-function and eT(k) = 

t2i 2 

+ Ur is the single particle energy and m*, Ur repre¬ 
sent the effective mass and the single particle potential, 
respectively, while k is the moment of the hole and q the 
external momentum transferred by the probe we use to 
excite the system. The latter can be taken along the z- 
axis without loss of generality [1^ . For simplicity, we will 
illustrate the case of SNM, but the formalism has been 
already generalized to the more general case of isospin 
asymmetric nuclear matter (ANM) [I^ |4^. Since the 
two Fermi surfaces are equal in SNM, we can drop the 
T index. The correlated RPA propagator is obtained by 
solving the Bethe-Salpeter (BS) equations as 


G(“)(ki,g,w) = (5aa'GffF(ki,(?,a;) -b Gffi;’(ki, q, w) ^ J \ki,k 2 )G°‘' {k 2 , q,u }), 

a' 


( 12 ) 


where a = {S, M, I) is a shorthand notation which stands 
for the quantum number of the system: the total spin 
S and its projection M along the z-axis and the total 
isospin I. Thanks to the particular form of the residual 
interaction \ we can solve the BS equation analyt¬ 

ically. The linear response function can be then obtained 
by integrating over the hole momentum ki as 

x(“)(Q,c^) = nd(G(“)(ki,g,cc)), (13) 

where = 4 is the degeneracy of the system. The ex¬ 
plicit expressions of x“ {q, cc) are quite complicated and 
we refer the reader to Ref. for their explicit expres¬ 
sions. To detect a pole we can simply look for the zeros 
of the expression 

l/x(^)(uj = 0,q) = 0, (14) 

for given values of transferred momentum q and den¬ 
sity p. In Fig. [7l we show in the (p, g)-plane the posi¬ 
tion of the instabilities obtained solving Eq. da. Apart 
from the spinodal instab ility which is related to the gas- 
liquid phase transition |46l | and manifests itself in the 
low-density region in the (0,0,0) channel, all the other in¬ 
stabilities are unphysical. In particular, we observe that 
all the functionals apart from Sill, present an instabil¬ 
ity in the spin channel close to the saturation density. 
Such an instability indicates a phase transition between 
spin-unpolarized to spin-polarized symmetric matter. In 


the present article, we investigate how these ferromag¬ 
netic instabilities are related to the problems previously 
detected in finite nuclei. 

As already discussed in Ref. [i^, the position of the 
poles can be different when we vary the isospin asym¬ 
metry of the system. In Fig. [51 we show the evolu¬ 
tion of the poles in asymmetric matter in the S' = I 
channel and different asymmetries for the BSk27 func¬ 
tional. We observe that changing the asymmetry pa¬ 
rameter Y = from 0 (i.e. SNM) to 0.2 (the bulk 

asymmetry of ^°®Pb) the vector-isovector instability gets 
closer to saturation density, and one could thus expect 
the appearance of instabilities in very neutron rich nuclei. 
As one can see from Fig. |8l the variation of the poles is 
roughly monotonic from SNM to PNM (i.e. Y = I) (see 
also discussion in Ref. El)- We remind that in this case 
both the isoscalar and isovector coupling constants keep 
their nominal value. For such a reason, this result can 
not directly compared with the one obtained in FigjS] 
Anticipating the result of the next section, we conclude 
that the stability criterion derived here should be applied 
for both SNM and PNM as discussed in Ref. [^ . 

V. DETECTING INSTABILITIES 

In Fig. 0 we illustrate the process used to extract the 
critical density pc at which a pole occurs in SNM: we 
consider the lowest density value at which the pole occurs 
in a given spin-isospin channel. Although we are not 







p (fm ) p (fm ) p (fm 



FIG. 7: (Color online) Instabilities in SNM for the functionals considered in the present article. The dashed-dotted horizontal 
line stands for the saturation density of the functional. See text for details. 
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M=0 M=1 



TABLE III: Critical coupling constant and critical density 
for the functionals studied in the present article extracted at 
n —>■ oo. The error bars we have added come from the uncer¬ 
tainty in the extrapolation we used for studying the variation 
of the critical point as a function of the shell number. For 
each interaction we also give the position of the pole for the 
nominal value of the coupling constant. 



isoscalar 

isovector 


pcjpsat 

Ppole / psat 

pc/psat 

Ppole j Psat 

Sill 

1 1 Q + U-iS4 

-L.-LO_q Qg 

3.29 

i.UO_Q Qg 

2.66 

SLy5 

l-OltS'.fe 

1.07 

0 07 + 0-20 
'j-y'_o.o5 

2.07 

T44 

0.961°;^? 

0.52 

u.yo_o.o8 

0.88 

BSk27 

Ill +0.33 
-*-•-*-^--0.10 

1.12 

i.UU_0.04 

2.48 

SAMi 

0 Q7+0-21 

u.y/_o.o6 

1.06 

(M O 

d d 
+ 1 
o 
o 

2.80 


FIG. 8: (Color online) Evolution of the critical densities in 
the S=1 channel in asymmetric matter for different values of 
the asymmetry parameter Y for the BSk27 functional. 




0.6 0 . 


FIG. 10: (Color online) Critical density pc/psat- The un¬ 
certainty band comes from the numerical extraction and the 
number of shells employed in the calculations. Solid bands 
refer to the isoscalar channel, while dashed bands refer to the 
isovector one. The vertical black lines refer to the best value 
as given in Tab lllll 


FIG. 9: (Color online) Position of the poles in the (p, (jt)-plane. 
The minimum of this curve defines the point (pcQc). 


presently interested in knowing the associated value of 
the transferred momentum for a given p^, it is anyhow 
useful to distinguish between long-wavelength and finite- 
size instabilities. Since we consider functionals originated 
from zero-range Skyrme interactions, there is no upper 
limit to the value of the exchange momentum between 
ph pairs. 

In Tab. uni we report the values of the critical densities 


using the coupling constants extracted from Tab.|TTl Each 
value of Pc is obtained using the best value of the critical 
coupling constant extracted in TabHIl the error bars are 
extracted using the higher and lower value respectively 
for both isoscalar and isovector channels. 

This value is shown in Fig llOl by means of a vertical 
solid black line in both isoscalar and isovector channels. 
We observe that most of the error bars fall in a ran ge o f 
the density pc € [psatA-'^'Psat]- Compared to Ref. [I8l |. 
it seems that performing RPA calculations allow us to 
explore slightly higher density regions. 

We can conclude that the stability criterion derived in 
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Ref. [l^ is also valid for the vector part of the functional 
and thus can be used to detect systematically this kind 
of problems. We also stress that the present method can 
be easily incorporated into a fitting procedure to avoid 
by construction this kind of problems, but in this case it 
is better to consider safer criteria by excluding a larger 
density interval. 

VI. CONCLUSIONS 

We have investigated the appearance of finite-size in¬ 
stabilities in RPA calculations of finite nuclei. The in¬ 
stabilities manifest trough the appearance of imaginary 
phonons and by an anomalous accumulation of strength 
in the low-energy part of the vibrational spectrum. By 
using six representative Skyrme functionals, we have 
identified the terms of the functional which are a pos¬ 
sible source of these instabilities and given upper limits 
to the values of their coupling constants. We have also 
performed systematic RPA calculations in symmetric nu¬ 
clear matter and we have observed that the presence of 
such a kind of instabilities can be related to the appear¬ 
ance of spurious phase transitions in the infinite medium 
around saturation density. As done in Ref. [I^ , we have 
thus extracted a quantitative criterion to detect the pres¬ 
ence of finite-size instabilities. Since LR calculations in 


SNM and PNM are not time consuming, such a method 
can be easily included into a fitting protocol, thus avoid¬ 
ing such kind of problems by construction. An example 
has been already presented in Ref. [^ . 

Although the finite range interactions do not seem to 
manifest such kind of pathologies since the finite range 
works as an effective cut-off in momentum space. It 
would be probably interesting in the future to do a similar 
analysis especially for the functionals including explicit 
tensor terms [471 - 14^ . which have been source of several 
instabilities in the Skyrme case [50j |. 
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